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We investigate a class of reaction processes in which particles move ballistically and react upon 
colliding. We show that correlations between velocities of colliding particles play a crucial role in the 
long time behavior. In the reaction-controlled limit when particles undergo mostly elastic collisions 
and therefore are always near equilibrium, the correlations are accounted analytically. For ballistic 
aggregation, for instance, the density decays as n ~ t~^ with £ = 2d/(d+3) in the reaction-controlled 
limit in d dimensions, in contrast with well-known mean-field prediction £ = 2d/(d + 2). 
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Ballistic-controlled reaction processes 0, IE IE El IE 
IE IE IE IE EE] exhibit rich atypical behaviors, e.g. the 
persistent dependence of decay exponents on the spatial 
dimension d implying absence of the upper critical di- 
mension; not surprisingly, ballistic-controlled processes 
proved very challenging to theoretical treatments. The 
key such process is ballistic aggregation pj in which par- 
ticles merge upon collisions so that mass and momentum 
are conserved (energy is necessarily lost). This model 
arises in various contents, e.g., it mimics the merging of 
coherent structures (like vortices or thermal plumes) and 
accumulation of cosmic dust into planetesimals . The 
one-dimensional (ID) version has also an interesting con- 
nection with dynamics of shocks representing solutions of 
the inviscid Burgers equation [Tj, [T]| ■ Ballistic aggrega- 



average momentum p and velocity V scale as 



tion was first investigated in a pioneering paper [lj by 
Carnevale, Pomeau, and Young who argued that basic 
physical quantities behave algebraically in the long time 
limit, e.g., the density decays as 



n(t) ~ t 



2d 
d + 2 



(1) 



in d dimensions. To understand this result, one can use 
a rate equation dn/dt = —n/r for the density. The 
mean time r between collisions related to the root mean 
squared (rms) velocity V, radius R, and density through 
nVrR ~ 1. Mass conservation implies that the aver- 
age mass is m ~ nT 1 . Therefore R ~ wr x l d and 



^ = - n 2 VR d-i = -n 1+1 ' d V. 
dt 



(2) 



The particle of mass m is formed from m original particles 
(we measure mass in units of the initial mass and velocity 
in units of the initial rms velocity). Assuming velocities 
of those original particles uncorrelated we find that the 
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p~m"'~, V — p/m^n'. (3) 

Plugging JSJ into J5J and solving for n(t) yields (JTJ. 

Surprisingly the prediction £ = 2d/(d + 2) for the de- 
cay exponent — perhaps the most known result in the 
field of ballistic-controlled processes — is erroneous. It 
turns out that the mean-field assumption that velocities 
of original particles contained within a typical aggregate 
particle are uncorrelated is incorrect in any finite dimen- 
sion — only when d — > oo and velocities are orthogonal 
to each other with probability one, they are indeed un- 
correlated. The failure of the mean-field no-correlation 
assumption Jjyl has not been appreciated because the re- 
sulting formula £d = 2d/ (d + 2) is correct both for d = 1 
and d — oo. (No trivial explanation of the former as- 
sertion is known yet the relation to the Burgers equa- 
tion via the particles<S=>shocks mapping 0, ^3| and the 
^2/3 g row th. of the separations between adjacent shocks 
established by Burgers many years ago |12| prove that 
£i = 2/3). Since £d monotonously increases with d, it 
is not surprising that the actual values are not so dif- 
ferent from the mean-field prediction QJ. Therefore the 
observed disagreement in two dimensions 7] could be at- 
tributed to insufficient scale of the simulations. Inter- 
estingly, the beauty of ballistic aggregation in ID where 
the model admits an exact solution [E and exhibits a 
deep connection to the Burgers equation has supported 
the incorrect prediction Q in higher dimensions. 

The purpose of this article is twofold. First, we clar- 
ify the role of velocity correlations in the general case, 
where they lead to significant deviations from mean-field 
predictions. Second, we propose a procedure that allows 
an analytical treatment of correlations for virtually any 
ballistic reaction process in the reaction-controlled limit; 
in particular, this method gives exact decay exponents. 

The no-correlation assumption is generally wrong for 
all ballistic-controlled processes, so we first demonstrate 
this assertion for one particularly simple process. We 
choose a toy ballistic aggregation model in which all par- 
ticles are identical and when two particles moving with 



2 



velocities vi and V2 collide, they form an aggregate par- 
ticle moving with velocity v = Vj + v 2 . Compared to the 
original ballistic aggregation model, the toy model has 
a number of advantageous properties. First, the volume 
fraction decays indefinitely thereby driving the system 
into the dilute limit and justifying ignoring multiple colli- 
sions. Second, the mean free path n _1 grows faster than 
the inter-particle distance n~ x l d for d > 1. These two 
features indicate that for d > 1 the Boltzmann equation 
approach is exact at large times. 

For the toy model, 10 becomes dn/dt = —n 2 V and the 
supposed absence of correlations gives V ~ n- 1 ' 2 . Thus 
the mean-field argument implies n ~ with £ = 2 in- 
dependently on dimension d. Numerically we find that 
this universality does not hold: £ increases with dimen- 
sion and approaches the mean-field prediction only when 
d — y co. For instance, we find (with an accuracy better 
than 1%) 



1.33 when d = 1, 
1.55 when d = 2, 
1.65 when d = 3. 



(4) 



These results were obtained by solving the Boltzmann 
equation describing the toy model 



dP(v, t) 
dt 



=J dudw P(u, t) P(w, t)\u — w S(u + w — v) 



2P(v,t) / dwP(w,t)|v - w| 



(5) 



We have solved this equation numerically implementing 
a Direct Monte Carlo (DMC) simulation scheme (see e.g. 
(l4j for the general method, and 9] for an application 
to a ballistic-controlled reaction process). The idea is to 
rephrase Eq. (JSJ as a stochastic process. In each step two 
particles, say with velocities u and w, are selected at ran- 
dom among a population of N particles, and the reaction 
happens with a probability proportional to |u — w|. If 
the reaction has been accepted, a new particle of velocity 
u + w replaces two original particles, so the number of 
particles changes to N — 1. The time is incremented by 
(N 2 \u — wl)^ 1 , and the process is iterated again. This 
numerical scheme allows us to treat systems with initial 
number of particles of the order of 10 7 . The master equa- 
tion associated to this Markov chain is precisely (J5J so 
that we obtain the numerically exact solution of our prob- 
lem. The exponent values J3J significantly differ from the 
mean- field prediction £ = 2, and leave no doubt that the 
no-correlation assumption is wrong. 

In the Boltzmann equation l|5|l. the relative velocity 
|v — w| gives the rate of collisions and its non- linear char- 
acter makes analytical progress hardly possible. An old 
trick to overcome this difficulty is to replace the actual 
relative velocity by the rms velocity ^5|. This results in 
the Maxwell model that played an important role in the 
development of kinetic theory 0,0]. For the toy model, 
we have (hereafter the dependence on time is suppressed 



for ease of notation) 
1 <9P(v) 



V dt 



J duP(u) P(v - u) - 2nP(v) 



(6) 



Integrating JJjJ we find that the density n = J dw P(w) 



satisfies dn/dt 



2 V while nV 2 — J dw w 2 P(w) re- 



mains constant. Hence V = n~ x l 2 and £ = 2 showing 
that the mean-field no-correlation approach is essentially 
the Maxwell model in context of ballistic processes [18j . 
The Maxwell model is an uncontrolled approximation to 
the Boltzmann equation for the hard sphere gas and, not 
surprisingly, the exponents found within this approach 
are generally erroneous (see 19] for an alternative sim- 
plification, the so-called very hard particle approach). Of 
course, one could anticipate that the exponent £ = 2 
characterizes the Maxwell model without computations 
- the essence of the Maxwell model, that is the fact 
that collisions are completely random, assures that the 
no-correlation condition does hold. 

We now present an argument that emphasizes the role 
and importance of correlations between velocities of col- 
liding particles and applies to all ballistic-controlled re- 
action processes. The key point is to supplement an 
evolution equation for the mass density by an evolu- 
tion equation for the density of kinetic energy. For 
an arbitrary ballistic-controlled reaction process we de- 
note P{m, v, t) the joint mass-velocity distribution func- 
tion, and e = mv 2 the kinetic energy of a given parti- 
cle (for the toy model, we set m = 1). The evolution 
equations for the density n and kinetic energy density 
nE = J mv 2 P(m, v, t) dm dv = n(mv 2 ) read 



dn 
~dt 



d{nE) 
dt 



i(Ae 



(7) 



The mean free 
' 1 gives n x l d V ~ t~ x 



The first equation is just the definition of the time de- 
pendent collision frequency 1/r, the second additionally 
contains the kinetic energy (Ae) co n lost on average in 
a binary collision. In the scaling regime the quantities 
(Ae) co ii and E = {mv 2 ) exhibit the same time depen- 
dence, so the dissipation parameter a = (Ae) co n/P is 
asymptotically time independent. From Eqs. Q we get 
d\n(nE)/d\nn — a, or V 2 — nE 
path argument r _1 ~ nVR^ 1 ~ t 
for ballistic aggregation. Combining these two relations 
and the definition of £ we obtain £ = (l/d+a/2)^ 1 . Sim- 
ilarly for the toy model nV 2 ~ n a and nV ~ t~ Y leading 
to £ = 2/(1 + a). 

To use this formalism, we must precisely define the 
collisional average involved in (JJJ . An av erag e change of 
a quantity A(l, 2) in a binary collision is |21j 

/AM _ Mrf2| Vl -v 2 r[A^(l,2)]P(l)P(2) 
^ coll ~ /dld2|vi-V2|"P(l)P(2) ' 1 ' 

where we have used shorthand notations i — (m,, Vj) and 
di = dmi dwi (i = 1,2). In the key case of hard spheres 
we have v = 1, whereas the cases v = and v = 2 



3 



correspond the Maxwell and very hard particle models 
[T^|. respectively. We now illustrate the formalism for 
the toy model. The kinetic energy lost in a collision is 
Ae = v 2 + v| — (vi + v 2 ) 2 = — 2 Vi ■ v 2 . Hence 



(Ae) 



/ dvi dv 2 |vi - V2r(vi ■ V2) f (vi) f (V2) 

/ dVi dv 2 |V! -V2| K P(vi)P(v 2 ) 



For the Maxwell model (1/ = 0), the isotropy of P(v, i) 
shows that (Ae) co ji = 0, so a = and £ = 2/(1 + a) — 2 
in agreement with our previous calculation. Similarly for 
very hard particles (y = 2, see 0) we use isotropy to 
simplify a and arrive at 



a 



J dvi dv 2 (vi ■v 2 ) 2 P(vi)P(v2) 
[Jrfv^P(v)] 2 



The isotropy allows to compute the ratio of the integrals 
to yield a — 2/d leading to £ = 2d/(d + 2). For other 
values of v, including the case of interest v = 1, the dis- 
sipation parameter a depends on details of the velocity 
distribution, and isotropy alone is not sufficient to de- 
termine a. The reason for the failure of the mean-field 
argument — which amounts to the complete neglect of 
collisional correlations ((vi • V2) C oll — 0) — is now clear: 
in general, a collision involving a pair (vi, v 2 ) with a neg- 
ative product Vi • v 2 < has a higher probability than a 
collision with Vi • v 2 > 0. The dissipation parameter a 
is therefore positive so that £ = 2/(1 + a) < 2. Thus the 
mean-field prediction £ = 2 is an upper bound for £. 
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FIG. 1: Density versus time in the 2D aggregation model. 
The continuous straight line has slope —0.86. The inset shows 
that the average energy {mv 2 } decays as t~ ' 28 . Mean-field 
predictions are shown by the dashed lines (slope —1 in the 
main graph and in the inset). The non- linear Boltzmann 
equation describing ballistic aggregation has been solved by 
DMC for a system of N = 4 x 10 7 particles. The initial 
density is no and u>o denotes the initial collision frequency of 
the equilibrium hard sphere fluid. 



The above framework applies to any irreversible pro- 
cess with ballistic transport. For ballistic aggregation, 
the omission of collisional correlations amounts to setting 
a = 1, i.e. that the typical energy dissipated in a col- 
lision is the mean kinetic energy per particle. However, 
particles with larger velocities undergo more frequent col- 
lisions so that the mean energy dissipated exceeds the en- 
ergy of a typical particle. Hence a — (Ae) co n/E > 1 so 
that £ is smaller than the mean-field prediction 2d/(d+2). 
Previous Molecular Dynamics (MD) simulations have 
shown that £ ~ 0.85±0.04 in 2D for low volume frac- 
tions, with scaling laws extending over 2 decades in time 
■ The DMC technique allows to reach much larger time 
scales. Figure H shows that after an initial transient, the 
density exhibits a clear power law behavior over 5 decades 
in time. We estimate £ ~ 0.86 ± 0.005, in agreement 
with MD simulations. The inset displays the behavior of 
E = {mv 2 ), the quantity that is (asymptotically) time 
independent according to the mean-field prediction yj; 
we find E ~ t" 28 We have also performed DMC 

and MD simulations in 3D giving £ ~ 1.06 ±0.01. As ex- 
pected, the actual values of £ are smaller than the mean- 
field prediction f = 2d/(d + 2)[23|. 

Ballistic-controlled processes are generally intractable 
analytically. Following the fruitful line of attack on dif- 
ficult problems — generalize them! — let us consider a 
process in which colliding particles react with probabil- 
ity e and scatter elastically with complementary proba- 
bility 1 — e. The mean-field no-correlation argument is so 
general that it applies to these processes; in particular, 
according to mean-field the exponent £ is independent 
on e. Remarkably, we can now compute the exponent 
£ for one special value of e, viz. for e — > + . In this 
reaction- controlled limit particles undergo mostly elastic 
collisions. Therefore, the particles are always at equilib- 
rium, i.e., the velocity distribution is Maxwellian. This 
key feature makes the problem tractable. Consider for 
instance the toy model. One can compute 



(Ae) 



/ dvi rfv 2 |vi - v 2 |"(vi ■ v 2 ) P(vi) P(v 2 ) 
/ dvi rfv 2 I vi - v 2 |^P( Vl )P(v 2 ) 



for arbitrary v when P(v) is Maxwellian |24|. In partic- 
ular for the important case v = 1 we obtain a = 1/d, so 
that £ = 2d/(d+ 1). This exact result provides a useful 
check of numerical scheme (our DMC simulations in two 
and three dimensions are indeed in excellent agreement 
with the theoretical prediction). In contrast, the mean- 
field no-correlation argument predicts £ = 2 irrespective 
of the value of e. We see again that this is correct only in 
the d — > 00 limit. Note also that £ = 2dj (d + 1) which is 
exact for the reaction-controlled (e — > + ) version of the 
toy model provides a better 'guess' for £ in the original 
(e — 1) model than the mean-field approach [compare 
£ = 1, 4/3 and 3/2 in ID, 2D, and 3D to the numerical 
values (3J|]. 

Remarkably, the exponent £ in the reaction-controlled 
limit of ballistic aggregation can be computed even 
though the mass distribution n(m) = n _1 J dv P(m, v) 
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is unknown. The important point is that when e — > + , 
the joint mass/kinetic energy distribution function fac- 
torizes. Then one finds a = 1 + 1/d, or equivalently 
£ = 2d/(d + 3) independently on II(m) |25j. This exact 
result of course agrees with DMC simulations. Interest- 
ingly, it also provides a reasonable approximation of £ for 
the original (e = 1) aggregation model: £ = 0.8 in 2D and 
1 in 3D, to be compared to 0.86 and 1.06, respectively. 

Many other ballistic-reaction processes are solvable in 
the reaction-controlled limit. For instance for ballistic 
annihilation 3] , there is no exact solution in any dimen- 
sion yet in reaction-controlled limit, the exact value of 
the density decay exponent is given by £ = 4d/(4d + 1). 
This result is in surprisingly good agreement with nu- 
merical values for e = 1: £ = 4/5 vs. 0.804 26] in ID; 
£ = 8/9 ~ 0.89 vs. 0.87 |3 in 2D; £ = 12/13 ~ 0.92 
against 0.91 in 3D. We have studied several other 
ballistic-reaction processes |20j, e.g., a simplified ballis- 
tic aggregation model in which mass and momentum are 
conserved yet the radius does not grow. For this model, 



the mean-field prediction is £ = 2/3 independently on 
dimension d, whereas in the reaction-controlled limit, we 
get the exact result £ = 2d/{2>d + 1) (i.e. 0.571 in 2D 
and 0.6 in 3D). It is again instructive to compare these 
values with numerical results for e = 1: £ ~ 0.60 in 2D, 
and £ ~ 0.62 in 3D. 

We have shown that correlations between velocities of 
colliding particles govern the behavior of all reacting pro- 
cesses with ballistic transport. We illustrated importance 
of correlations on several models and demonstrated that 
ignoring correlations is equivalent to using the Maxwell 
model, which is an uncontrolled approximation of the 
hard-sphere gas. We also devised a procedure that clari- 
fies the role of correlations in the general case and allows 
an exact computation of decay exponents in the reaction- 
controlled limit, when particles undergo mostly elastic 
collisions and therefore are always near equilibrium. The 
failure of mean-field theory to describe this limit empha- 
sizes inevitable presence of correlations in all reacting 
processes with ballistic transport. 
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